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Precision measurements of the anomalous magnetic moment of the muon depend on the correct 
collective mode assignment of the coherent betatron oscillations in the beam. Presently, there is a 
disagreement between experiment and theory for the horizontal betatron frequency. The discrepancy 
is here resolved by our computations of electrostatic beam focusing. The correct treatment of the 
betatron effects renders less likely the need for non-standard model corrections to the muon g-2. 



PACS numbers: 29.20.Fj, 29.27.Bd, 29.27.Fh, 13.35.Bv, 13.30.Em 
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I. INTRODUCTION 



There has been considerable interest in precision 
measurements |l], @, [| |], ||, 0, frl of the muon anomalous 
magnetic moment factor 
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The anomaly factor is measured from the chiral asymme- 
try rotation frequency 
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Largely due to experimental improvements in creating 
an extremely homogeneous magnetic field B, there have 
been recent improvements Q in the accuracy of k deter- 
minations. 

In addition to the anomaly frequency of Eq.(2), there 
are three other frequency scales of interest in the exper- 
iments. The cyclotron frequency is given by 
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which is also determined by the homogeneous magnetic 
field. Due to electrostatic focusing capacitors, there are 
two more frequencies having to do with collective oscilla- 
tions of the beam as a whole. These are the vertical beta- 
tron oscillations at frequency f2ii and the horizontal beta- 
tron oscillations at frequency il±. In the literature^, |), 
the ratios (fiii/a; c ) and (£l±/u c ) are called, respectively, 
the vertical tune and horizontal tune. 

Recent experiments fji p| yield a vertical tune of 



with a measured effective index of 
n « 0.137 . 



(4) 



(5) 



Previous theoretical considerations 61 [9) lead to a hori- 
zontal tune of 



0.929 



(6) 



The horizontal tune has been measured by the muon de- 
cay oscillations themselves. (See and || wherein the 
notation k>& is used for f2j_.) 



0.0702 . 
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experiment 



The discrepancy between the previous theoretical Eq.(6) 
and the experimental Eq.(7) is evident. 

One purpose of our work is to show that the previous 
theory is at fault. The index n describes the electric field 
gradient inside the electrostatic focusing capacitors; It is 



IBI 



dE u 



dr\\ 



(8) 



where po ~ 711.2 cm is the radius of the cyclotron, 
v = (dr/dt) is the muon velocity, E\\ — (E-B/|B|) is 
the component of the electric field parallel to the mag- 
netic field and ri = (r • B/|B|). The averaging is over a 
cyclotron orbit. 

Some insight into the error of the previous theoretical 
Eq.(6) becomes evident from the fact that when the elec- 
tric field vanishes, lim Iw o(^-L/ a; c)expectcd = 1, the hori- 
zontal tune is unity. Such "unit horizontal tunes" are in 
reality not collective modes at all. A charged particle in 
a uniform magnetic field (and zero electric field) moves in 
a circle with a fixed center. If the center of circular orbit 
is not quite equal to the center of the cyclotron machine 
(which is always the case), then the above corresponds 
to what has been called a unit horizontal tune. 

However, the true nature of the horizontal betatron 
oscillation is that the center of the cyclotron orbit also 
rotates if and only if a non-zero electric field is applied. 
A central result of this work is that the center of the 
cyclotron orbit also moves in a circle at a horizontal tune 
of 



(5) 



0.0685 



(9) 



expected 



Our theoretical Eq.(9) is in satisfactory agreement with 
the experimental Eq.(7). 

In order to understand the nature of the horizontal 
betatron tune, consider a toy "hula hoop" popular with 
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children near the middle of the last century. Pretend 
there are particles rotating within (and relative to) the 
hoop with angular velocity u> c . Further imagine the hoop 
rotating as a whole "off center" with an angular velocity 
f2_i_. Now remove the hoop since it is imaginary anyway. 
The remaining particle motions constitute a beam with 
a horizontal tune. 

In Sec. II, the dynamics of the muon motion will be 
discussed starting from the Hamiltonian 

H = v/c 2 (p - (e/c)A) 2 + c A M 2 + e$, (10) 

where the magnetostatic and electrostatic fields are 
given, respectively, by 

B = curlA and E = grad$. (11) 

For simplicity of presentation, we shall employ quantum 
mechanics for describing the beam dynamics but we shall 
not employ (for now) the spin dynamics which must be 
treated with a Dirac Hamiltonian and requires somewhat 
long and delicate spinor algebra|Io| Q. The full quan- 
tum Dirac spinor algebra will be discussed elsewhere and 
should serve as a more accurate quantum replacement for 
the classical BMT equations |i~2]] . 

In Sec. Ill, the coherent betatron frequencies will be 
derived from the Hamiltonian Eq.(10). A physical de- 
scription of the coherent betatron oscillations will be pro- 
vided. In the concluding Sec. IV, the consequences for 
the (g — 2) analyses due to the corrected horizontal tune 
theory will be explored. 



The components of the cyclotron radius vector in Eq.(15) 
do not commute. From Eqs.(14)-(16) 

/ ihc\ 

\*> ( 17 ) 

The vector giving the position in the horizontal plane of 
the cyclotron orbit center 

R= (R x ,Ry,0) (18) 

is defined so that 

r = rn + rj_ and r± = p + R. (19) 

The components of the cyclotron orbit center position 
also do not commute; i.e. 

[Rx > Ry] = ~{^) ■ (20) 

We note en passant as an example of a noncommutative 
geometry |13| that the center coordinates R of the cy- 
clotron orbit live in the noncommutative geometric plane 
in virtue of Eq.(20). Finally, we have the vanishing com- 
mutators: 

[p,R] = 0, [p,z] = [R,z} = 0, (21) 

and 

[PiPz] = [R ( Pa] = 0. (22) 



II. QUANTUM BEAM DYNAMICS 

The Hamiltonian in Eq.(10) may be considered to be 
an operator by employing the canonical commutation re- 
lations 

\pi,rj] = -ihSi-j. (12) 

We also consider a coordinate system in the laboratory 
frame in which the uniform magnetic field vector is along 
the z = 7™ 1 1 axis; i.e. 

A = (A x ,Ay,0) and B = (0,Q,-B). (13) 

From Eqs.(ll)-(13) it follows that 

Px - -A x , Py - -AyJ = ( — — J . (14) 

The cyclotron radius vector in the horizontal plane 

P={Px,Py,0) (15) 

points from the cyclotron orbit center to the charged par- 
ticle position. In mathematical terms 

^(^) BX ( P ^ A ) < 16 > 




FIG. 1: The electrostatic focusing field pushes on a particle 
in a cyclotron orbit of radius vector p. The force causes the 
instantaneous center of the cyclotron orbit to rotate at an- 
gular velocity il± with a radius R. As schematically shown 
above, the resulting wobble in the particle orbit represents 
the horizontal coherent betatron oscillation. 
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The Hamiltonian of Eq.(lO) now reads 

H = Vc^pl + ( e B) 2 p 2 + c 4 M 2 + e<P(p + R, z). (23) 

To calculate the rate at which R changes with time, one 
may use 

R=l[iT,R]=r^[*(p + R,*),R]. (24) 

Thus, the center of the cyclotron orbit drifts with a 
velocity which depends only on the electromagnetic fields 
acting on the particle; i.e. 



V = R = c 



E x B 

— w 
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The method of deriving Eq.(25) using quantum mechan- 
ics is due to Schwingerj 
In the direction para 



14J. 

lei to the magnetic field, 
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where 



£ = ^c 2 P 2 z + (eB) 2 p 2 + c A M 2 = 1 Mc 2 . 
Furthermore, 



Pz = -[H,p z ] = ij) [$(p + R,*),p» 



Eqs.(27) and (28) imply 



M- 



dt 



= eE, 



(26) 



(27) 



(28) 
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If the Hamiltonian in Eqs.(lO) and (23) were treated 
using classical mechanics, then the operator equations 
of motion would turn into equations for particle orbits. 
There is sometimes simplicity in orbital pictures. Note 
that the classical orbital equations are recovered by re- 
placing commutators with Poisson brackets: 



(quantum) — [a, b] 
n 



(classical) {a, b} 



(30) 



where 

da \ /db\ _ /dy 
dv) \dv) V<9p, 

In the classical orbit view wherein 7 = {H, 7}, one finds 



{a,b} 



Or 



(31) 



.»/<•- [ ^ ) s i:. 



(32) 



The quantum length scale for which the classical orbit 
picture breaks down can be deduced from the noncom- 
mutative geometry of the cyclotron center coordinate R. 
From Eq.(20), the quantum uncertainties obey 



AR X AR V > (L 2 /2), 



(33) 



where 



L=y/(hc/\eB\ 



(34) 



For the muon experiments of interest L « 2.0 x 10 -6 cm, 
which allows for a classical orbit picture of high accuracy. 
Even in the classical case, the decomposition r± = p + H 
of the position in the plane is crucial for understanding 
the horizontal betatron oscillation mode. 



III. COHERENT BETATRON OSCILLATIONS 

To compute the betatron vertical and horizontal tunes, 
one may expand the electric field E((0 + R, z) in powers of 
the small displacements z and R. Inside the electrostatic 
focusing capacitors 



Defining 
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(36) 



and taking only first powers of the displacement z in 
Eq.(29) yields 



(37) 



M7 [ — ) z= -Qz. 



To this accuracy, the vertical betatron oscillation fre- 
quency is given by 



M 7 f2f| = Q 



(38) 



in agreement with Eqs.(4) and (8). 

Expanding the electric field in Eq.(25) to linear order 
in R yields 



dR 

~dt 



= (^)< R * B > 



(39) 



where Eq.(36) has been invoked. In the z = plane, the 
vector R rotates with angular velocity whose magnitude 
is given by 



cQ 



2eB 



(40) 



in agreement with Eqs.(8) and (9); i.e. from Eqs.(3), (38) 
and (40) it follows that 



(41) 



The horizontal oscillation is depicted in FIG. 1. If 
a particle is injected into the cyclotron at a somewhat 
skewed angle, then the center of the cyclotron orbit will 
be somewhat removed from the center of the cyclotron 
machine. The instantaneous center of the cyclotron orbit 
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will itself move in a small circle at an angular velocity Q±. 
The small circular rotations of R induce a wobble in the 
particle's (almost) circular orbit, at the same frequency 
f2j_. If a particle is injected into the cyclotron with a 
small velocity component v z normal to the plane, then 
the height z of the particle undergoes oscillations with 
frequency fin. Since, in virtue of Eq.(36), there is but 
one electric field gradient parameter Q, the two betatron 
frequencies are not independent. They are related by 
Eq.(41). 



IV. CONCLUSION 

On a phenomenological basis, the chiral asymmetry of 
the muon decay is fit to experimental data by a func- 
tion of the product form 

Ptot (t) = P K (t)P h 

etatron (t)Pioss{t), (42) 

where the chiral rotation function is parameterized by 

P*(t) = e~ t/7r " (1 + A K smfat + cf> K )) , (43) 

the coherent horizontal betatron oscillation is parameter- 
ized by 

Pbetatron{t) = ^ (1 + A b COs(fi ± i + (j) b )) , (44) 

and the extra loss of muons from unknown physical pro- 
cesses are parameterized by 

Plo SS {t) = l + nje-*/*!. (45) 

The agreement between the experimental best fitB 



2.05, 



(46) 



n / experiment 

and our theoretical prediction based on Eqs.(8) and (9) 



2.01, 



(47) 



gives us confidence that we have properly identified the 
source of horizontal betatron oscillations. It is something 
of a mystery to us why the experiment was designed so 
that fl± is so very close to being a "second harmonic" of 
the anomaly frequency u> K . 

If one wishes to go beyond the phenomenological view 
of how the electromagnetic field F al , — d u A v — d v A u 



affects the coupling g = 2(1 
Dirac-Schwinger equation 



n), one must employ the 



-iH^du + Mc ■ 



( I" ,; \ T/ „. , 
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a^F^u >fl> = 0, (48) 



where = (d u — (ie/hc)A u ). The interaction between 
the electromagnetic field and the anomalous magnetic 
moment is thereby 



W/tuF"" = k(B • a - i 75 E • a) 



(49) 



K / theory 



The electric field coupling of Eq.(49) via E(p + R, z) 
contains those degrees of freedom (p z , z, p x , p y ) which en- 
ter into the kinetic energy £ = Mc 2 ^ in Eq.(27) and also 
the cyclotron center coordinates (R x , Ry) which rotate as 
in Eq.(39). 

The experimentally observed coupling of the chiral ro- 
tation frequency lu k to the horizontal betatron frequency 
Qj_ is mainly due to the variation in R. This coupling to 
R dominates the electric field corrections to oj k due to p. 

In assessing the electric field corrections to w K Q], due 
to the variations in the radius (r ^ = p + R) of the beam 
orbit, the effect was doubly counted: (i) The time varia- 
tions in r_L were first counted in the function Pbetatron(t) ■ 
(ii) The time variations in rj_ were counted yet again as 
a BMT electric field correction. Leaving out the over 
counting of the electric field corrections makes the ex- 
perimental need for non-standard model corrections to 
the muon {g — 2) less compelling. 
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